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Abstract 

We consider the spin-i Ising chain in a regularly alternating transverse field to examine the effects of regular 
alternation on the quantum phase transition inherent in the quantum Ising chain. The number of quantum phase 
transition points strongly depends on the specific set of the Hamiltonian parameters but never exceeds 2p where p 
is the period of alternation. Calculating the spin correlation functions numerically (for long chains of up to 5400 
sites) and determining the critical exponents we have demonstrated that two types of critical behavior are possible. 
In most cases the square-lattice Ising model universality class occurs, however, a weaker singularity may also take 
place. 
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We consider the spin-^ Ising chain in a transverse 
field [1,2,3] assuming that the Hamiltonian parameters 
(i.e. the nearest neighbor interactions and the on-site 
fields) vary regularly along the chain with a finite pe- 
riod of alternation p. The Hamiltonian of the model 
reads 

H = ^2 2I n s x n s* l+1 + ^2 n n s z n (1) 

n n 

and the sequence of parameters in (1) is 

7if2l/2^2 • ■ ■ Ip^lpllflll2^2 ■ ■ ■ IpQp .... 

Recently, it has been shown [4] that a regularly alter- 
nating transverse field may strongly influence the quan- 
tum phase transition inherent in the uniform chain 
[1,2,3]. In what follows we consider the chain of period 
2 with fii^ = H± An and J„ = — 1. Then the quantum 
phase transition points follow from the condition [5] 



(n* + An) (n* - An) = ±1. 



(2) 
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= zt VAn 2 + i if 
= ±VAn 2 ± i if 



Eq. (2) gives two critical fields n* 
An < 1 or four critical fields n* 
An > 1. Moreover, if An = 1 in addition to n* = 
we have one more critical field n* = 0. Further 
we restrict ourselves to the case An = 1 and examine 
the critical behavior in the vicinity of the two repre- 
sentative critical points n* = \fl and n* = 0. This 
is by no means a restrictive example since only two 
types of critical behavior are expected for the model 
(1) [4] . The critical behavior is characterized by a set of 
critical exponents. Namely, the order parameter (Ising 
magnetization) vanishes as m x ~ (n* - The two- 
point spin correlation function (s£s£ + 2 r ) — (s%) (s^ + 2 r } 
in the limit r — > oo decays as ~ (2r)~ v exp 2pV 

where £ ~ |Q — Q*\~ v . The transverse static suscepti- 
bility diverges as x z ~ |n — n*|~ Q and the energy gap 
disappears as A ~ n — n* \ vz . From the analytical cal- 
culations [4] we know that in the vicinity of n* = \/2 
we have vz = 1 and a = (x z exhibits a logarithmic 
singularity) whereas in the vicinity of n* = we have 
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Fig. 1. Power-law decay of spin correlations at the critical 

points. Main plot: (Sn s n+2r) vs (log-log plot, N — 2000, 

n — 920, O* — (up), fl* — y/2 (down), the solid lines are 

l _ l - 

0.16 (2r) * and 0.08 (2r) Inset: vs (Q = (down), 

£1* = V2 (up), the solid curves arc 0.425iV~i and 0.467V~H). 
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Fig. 2. The order parameter behavior in the critical regions 

(N = 2000: empty symbols, N = 5400: full symbols) [4]. The 

i / x I 

solid curves arc 0.5f23 (a) and 0.5 (V5 - Q) 8 (b). 

^2 = 2 and a — — 2 (%^- exhibits a logarithmic singu- 
larity) . In the present paper we complete these analyti- 
cal findings by precise numeric results for other critical 
exponents rj, (3 and v computing two-point spin corre- 
lation functions for chains of a few thousand sites [6]. 

We start our analysis considering {s^ l sf l+2r ) at the 
critical fields fi* = and Q,* — \/2. The results re- 
ported in the main plot in Fig. 1 clearly indicate the 
power-law decay of spin correlations with r\ = \ at both 
critical points Q.* = (squares) and fi* = \[2 (dia- 
monds). We calculate the on-site x-magnetization mj, 
mj 2 = lim r ^ 00 (sjsj +2r ), for finite chains of N sites 
putting j = \N, 2r = |iV. mj vanishes at both criti- 
cal fields as as can be seen in the inset in Fig. 1 
(squares refer to fi* = 0, diamonds refer to Q.* — \/2). 

We turn to the critical region. The results for m x — 
| (mf + mf) are shown in Fig. 2 (see also [4]). In the 
vicinity of fi* = (Fig. 2a) we observe a power-law 
vanishing of the order parameter with f) = \. In the 




0.1 1.4 1.41 1.42 1.43 n 



Fig. 3. The inverse correlation length behavior in the critical 
regions (TV — 2000: empty symbols, N — 5400: full symbols). 
The solid curves arc U 2 (a), 3\Sl-V2\, n < V2 and V^|fJ-v^|, 

n > V2 (b). 

vicinity of Q* = \/2 (Fig. 2b) we observe a power- 
law vanishing of the order parameter with (3 — | . Fi- 
nally, in Fig. 3 we present the results for the correlation 
length. While in the vicinity of Q* — we find v = 2 
(Fig. 3a) in the vicinity of fi* = \/2 we clearly observe 
v = 1 (Fig. 3b). 

To summarize, combining the analytical and numer- 
ical results we have found that the spin- \ Ising chain in 
a regularly alternating transverse field may exhibit the 
critical behavior which belongs to the square-lattice 
Ising model universality class with ( 3=|,z/=l, 7) = 
j, a — 0, z — 1 (as the considered chain at Q* — \/2) 
and a weaker singularity with (3 = \, v = 2, 77 = i, 
a — —2, z = 1 (as the considered chain at Q* = 0). 
Let us finally note, that the spin correlation functions 
required for the estimation of exponents j3, v, rj should 
be possible to calculate analytically as in the uniform 
case. We left this problem for future studies. 
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